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Abstract 

We derive the non-universal gaugino mass ratios in a supergravity (SUGRA) frame- 
work where the Higgs superfields belong to the non-singlet representations 54 and 770 
in an S'O(IO) Grand Unified Theory (GUT). We evaluate the ratios for the phenomeno- 
logically viable intermediate breaking chain SU{A)c x SU(2)l x SU(2)r (G^)- After 
a full calculation of the gaugino mass ratios, noting some errors in the earlier calcula- 
tion for 54, we obtain, using the renormalisation group equations (RGE), interesting 
low scale phenomenology of such breaking patterns. Here, we assume the breaking 
of the SO(W) GUT group to the intermediate gauge group and that to the Standard 
Model (SM) take place at the GUT scale itself. We also study the collider signatures in 
multilepton channels at the Large Hadron Collider (LHC) for some selected benchmark 
points allowed by the cold dark matter relic density constraint provided by the WMAP 
data and compare these results with the minimal supergravity (mSUGRA) framework 
with similar gluino masses indicating their distinguishability in this regard. 
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1 Introduction 



With the Large Hadron Collider (LHC) about to be operative shortly, the search for physics 
beyond the standard model (SM) has reached a new height of excitement. Low energy (TeV 
scale) supersymmetry (SUSY) has persistently remained one of the leading candidates among 
scenarios beyond the Standard Model (SM), not only because of its attractive theoretical 
framework, but also for the variety of phenomenological implications it offers [1, 2, 3, 4]. 
The stabilization of the electroweak symmetry breaking (EWSB) scale and the possibility of 
having a cold dark-matter (CDM) candidate with conserved .R-parity (R = (— i)( 3B + L + 25 ')) 
[5] are a few elegant phenomenological features of SUSY. Side by side, the possibility of 
paving the path towards a Grand Unified Theory (GUT) is one of its most exciting theo- 
retical prospects, where one can relate the SU(3), SU(2) and U(l) gauge couplings and the 
corresponding gaugino masses at a high scale [4, 6, 7]. 

The most popular framework of SUSY-breaking is the minimal supergravity (mSUGRA) 
scheme, where SUSY is broken in the 'hidden sector' via gravity mediation and as a result, 
one can parametrise all the SUSY breaking terms by a universal gaugino mass (M 1 / 2 ), a 
universal scalar mass (mo), a universal trilinear coupling parameter Aq, the ratio of the 
vacuum expectation values (vev) of the two Higgses (tan/3) and the sign of the SUSY- 
conserving Higgs mass parameter, (sgn(fi)) [4, 8]. 

However, within the ambit of a SUGRA-inspired GUT scenario itself, one might find some 
deviations from the above-mentioned simplified and idealized situations mentioned above. 
For instance, the gaugino mass parameter (Mi/ 2 ) or the common scalar mass parameter 
(m ) can become non-universal at the GUT scale. In this particular work, we adhere to 
a situation with non-universal gaugino masses in a supersymmetric scenario embedded in 
SO(10) GUT group. 

Gaugino masses, arising after GUT-breaking and SUSY-breaking at a high scale, crucially 
depend on the gauge kinetic function, as discussed in the next section. One achieves universal 
gaugino masses if the hidden sector fields (Higgses, in particular), involved in GUT-breaking, 
are singlets under the underlying GUT group. However if we include the higher dimensional 
terms (dimension five, in particular) in the non-trivial expansion of the gauge-kinetic func- 
tion, the Higgses belonging to the symmetric products of the adjoint representation of the 
underlying GUT group can be non-singlets. If these non-singlet Higgses are responsible for 
GUT breaking, the gaugino masses Mi, M 2 and M 3 become non-universal at the high scale 
itself. It is also possible to have more than one non-singlet representations involved in GUT 
breaking, in which case the non-universality arises from a linear combination of the effects 
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mentioned above. 

Although this issue has been explored in earlier works, particularly in context of the 
SU(5) [9, 10, 11], there had been one known effort [12] to study the same in context of 
the SO(10). In this paper, we calculate the non-universal gaugino mass ratios for the non- 
singlet representations 54 and 770, based on the results obtained in [13], for the intermediate 
gauge group, namely, Pati-Salam gauge group SU(4) c x SU{2) L x SU{2) R (G422) with con- 
served .D-parity [14]. In order to understand the low-energy phenomenology of such high 
scale breaking patterns, which is indeed essential in context of the LHC and dark mat- 
ter searches, we scan a wide region of parameter space using the renormalisation group 
equations (RGE) and discuss the consequences in terms of various low-energy constraints. 
For example, we discuss the consistency of such low-energy spectra with the radiative elec- 
troweak symmetry breaking (REWSB), Landau pole, tachyonic masses etc. We also point 
out the constraints from stau-LSP (adhering to a situation with conserved i?-parity and 
hence lightest neutralino- LSP), and flavour constraints like b — > 37 for all possible com- 
bination of the parameter space points. Here, we assume the breaking of the 5*0(10) to 
the intermediate gauge group and that to SM takes place at the GUT scale. To study the 
collider signatures in context of the LHC, we choose some benchmark points (BP) consistent 
with the cold dark matter [5] relic density constraint obtained from the WMAP data [15]. 
We perform the so-called 'multilepton channel analysis' [16, 17] in same — sign dilepton, 
opposite — sign dilepton, trilepton, inclusive four — lepton channels associated with jets, 
as well as in hadronically quiet trilepton channel at these benchmark points. We compare 
our results with WMAP allowed points in mSUGRA, tuned at the same gluino masses 1 . 

There has been a lot of effort in discussing various phenomenological aspects [18, 19, 
20, 21, 22] of such high scale non-universality and its effect in terms of collider signatures 
[23, 24, 25], our analysis is remarkable in the following aspects: 

• Apart from noting some errors in the earlier available calculations of non-universal 
gaugino mass ratio for the representation 54, we present the hitherto unknown ratio 
for the representation 770 for the intermediate breaking chain £422- 

• While we discuss the consistency of the low-energy spectra obtained from such high 
scale non-universality in a wide region of the parameter space, we study the collider 
aspects as well in some selected BPs in context of the LHC. 

1 We are benchmarking a potentially faking SUGRA scenario in terms of the gluino mass, which has a 
very important role in the final state event rates. The corresponding sfermion masses are fixed by requiring 
the fulfillment of WMAP constraints. 
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• In order to distinguish such non-universal schemes from the universal one, we compare 
our results at these chosen BPs with WMAP allowed points in mSUGRA tuned at 
the same gluino masses. We identify a remarkable distinction between the two, which 
might be important in pointing out the departure in 'signature-space' [26, 27, 28, 29] 
in context of the LHC for different input schemes at the GUT scale. 

The paper is organized as follows. In the following section, we calculate the non-universal 
gaugino mass ratios. The low-energy spectra, their constancy with various constraints and 
subsequently the choices of the BPs have been discussed in section 3. Section 4 contains 
the strategy of the collider simulation and the numerical results obtained. We conclude in 
section 5. 

2 Non-universal Gaugino mass ratios for 50(10) 

Here we calculate the non-universal gaugino mass ratios for non-singlet Higgses belonging 
to the representations 54 and 770 under SO (10) SUSY-GUT scenario. 

We adhere to a situation where all soft SUSY breaking effects arise via hidden sector in- 
teractions in an underlying supergravity (SUGRA) framework, specifically, in SO (10) gauge 
theories with an arbitrary chiral matter superfield content coupled to N=l supergravity. 

All gauge and matter terms including gaugino masses in the N=l supergravity Lagrangian 
depend crucially on two fundamental functions of chiral superfields [30]: (i) gauge kinetic 
function f a p(<&), which is an analytic function of the left-chiral superfields $j and transforms 
as a symmetric product of the adjoint representation of the underlying gauge group (a, /3 
being the gauge indices, run from 1 to 45 for 5*0(10) gauge group); and (ii) G(<&i, $*), a real 
function of $j and gauge singlet, with G = K + ln\W\ (K is the Kahler potential and W is 
the superpotential) . 

The part of the N=l supergravity Lagrangian containing kinetic energy and mass terms 
for gauginos and gauge bosons (including only terms containing the real part of /(<&)) reads 

e- l C = ~Ref^(4>)(-l/2\ a p\^) - htef^F^F^ 

+\e~ G/2 G\(G- 1 )mi^)/d^}^ + h.c (1) 

where G l = dG/d(pi and (G -1 )* is the inverse matrix of G J j = dG / d<p* l d(f)j, \ a is the 
gaugino field, and <j) is the scalar component of the chiral superfield and F" is defined in 
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unbroken SO (10). The F-component of <£> enters the last term to generate gaugino masses 
with a consistent SUSY breaking with non-zero vev of the chosen F, where 

F* = ^e-^miG-'m (2) 

The & s can be a set of GUT singlet supermultiplets $ 5 , which are part of the hidden 
sector, or a set of non-singlet ones $ N , fields associated with the spontaneous breakdown of 
the GUT group to SU(3) x SU (2) xU(l) . The non-trivial gauge kinetic function f a p(&) can 
be expanded in terms of the non-singlet components of the chiral superfields in the following 

way 



faeiv) = fo(^ s )s a , + £ + O(^) 2 (3) 

where fo and £ N are functions of chiral singlet superfields, essentially determining the 
strength of the interaction and M is the reduced Planck mass= M P i/\/8n. 

In Equation (3), the contribution to the gauge kinetic function from Q N has to come 
through symmetric products of the adjoint representation of the associated GUT group, 
since f a p on the left side of Equation (3) has such transformation property for the sake of 
gauge invariance. 

For SO (10), one can have contributions to f a p from all possible non-singlet irreducible rep- 
resentations to which <& N can belong : 



(45 x A5) symm = 1 + 54 + 210 + 770 (4) 

As an artifact of the expansion of the gauge kinetic function f a p mentioned in Equation (3), 
the gauge kinetic term (2nd term) in the Lagrangian (Equation (1)) can be recast in the 
following form 

ReMflF^F^ = ^p-Tr(F^ N Fn (5) 

where F^ V) under unbroken 5*0(10), contains U(l)y, SU(2) L and SU(3)c gauge fields. 

Next, the kinetic energy terms are restored to the canonical form by rescaling the gauge 
superfields, by defining 

F% v ^F« u = (Ref a ^F^ u (6) 
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and 

\° ^ \° = (Ref^hP (7) 
Simultaneously, the gauge couplings are also rescaled (as a result of Equation (3)): 

g a (M x ){Ref a ^S aP = g c {M x ) (8) 

where g c is the universal coupling constant at the GUT scale (M x )- This shows clearly that 
the first consequence of a non-trivial gauge kinetic function is non-universality of the gauge 
couplings g a at the GUT scale [9, 10, 31, 32]. 

Once SUSY is broken by non-zero vev's of the F components of hidden sector chiral 
superfields, the coefficient of the last term in Equation (1) is replaced by [9, 10] 

(F^) = 0(msM) (9) 

where ms = exp{—^-) is the gravitino mass. Taking into account the rescaling of the 
gaugino fields (as stated earlier in Equation (7)) in Equation (1), the gaugino mass matrix 
can be written down as [9, 11] 

M a {M x )5 a (3 = ^ — — - (10) 

which demonstrates that the gaugino masses are non-universal at the GUT scale. 

The underlying reason for this is the fact that {f a p) can be shown to acquire the form 
fa$af3i where the f a 's are purely group theoretic factors, as we will see. On the contrary, if 
symmetry breaking occurs via gauge singlet fields only, one has f a p = foS a /3 from Equation 
(3) and as a result, (f a p) = fo- Thus both gaugino masses and the gauge couplings are 
unified at the GUT scale (as can be seen from Equations (8) and (10)). 

As mentioned earlier, we would like to calculate here, the f a 's for Higgses (& N ) belonging 
to the representations 54 and 770 which break 5*0(10) to the intermediate gauge group 
SU(A) C x SU(2) L x SU(2) R with unbroken D-parity (usually denoted as G 42 2p) 2 ■ We 
associate the non-universal contributions to the gaugino mass ratios with the group theoretic 
coefficients f a 's that arise here. It in turn, indicates that we consider the non-universality 
in the gaugino masses of 0(1). This, however, is not a generic situation in such models. This 
can be rather achieved under some special conditions like dynamical generation of SUSY- 
breaking scale from the electroweak scale, no soft-breaking terms for the GUT or Planck 

2 Higgs, that breaks G422P to SM, doesn't contribute to gaugino masses. 
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scale particles and with the simplified assumption Mqut = Mpi which is also reflected in 
the RGE specifications. 

The representations of SO (10) [33], decomposed into that of the Pati-Salam gauge group 

are 

SO(10) ->■ G422 = SU(A) C x SU(2) L x SU(2) R 

45 ->■ (15, 1,1) ©(1,3, 1)©(1, 1,3) ©(6, 2, 2) 

10 (6, 1,1) ©(1,2, 2) (11) 

16 -+ (4,2,1) © (4,1,2) 

Using the 5*0(10) relation, (10 x 10) = 1+45 + 54, one can see that vev of 54-dimensional 
Higgs ((54)) can be expressed as a (10 x 10) diagonal traceless matrix. Thus the non-zero 
vev of 54-dimensional Higgs can be written as [32] 

< 54 >= diag(3, 3, 3, 3, , -2, -2, -2, -2, -2, -2). (12) 

2v 15 



Since (45 x 45) s?/m = 1 + 54 + 210 + 770, one can write the non-zero vev [13] of 770- 
dimensional Higgs as (45 x 45) diagonal matrix: 

(770) = diag(-A, , -4, -10, -10, 5, , 5) (13) 

V180 v ' ' v ' 

15 3+3 24 

In the intermediate scale (M c ), G 42 2 is broken to SM gauge group. Here SU(A) C is 
broken down to SU(3)c x U(1)b-l an d at the same time, SU(2)r is broken to U(1)t 3 r- It 
is noted that SU(2)r x SU(4)c is broken to SU(3)c x U(l)y and hence the hypercharge is 
given as Y = T 3R + ^(B — L). Below we note the branchings of SU(A)c representations: 

SU(A) C -> SU(3) C x U(1) B _ L 

4 ->■ 3 1/3 ©l_i 

15 ->■ 8 © 3 4/3 © 3_ 4 /3 © lo (14) 

10 ->■ 62/3 ©3_ 2 /3 © 1-2 

6 — > 3_2/3 © 32/3 

Combining these together, we achieve the branchings of 50(10) representations in terms 
of the SM gauge group. 

5O(10) : SU(3) C x SU(2) L x U(1) Y 
45 
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(15,1,1) 



(1,3,1 
(1,1,3 
(6,2,2 

10 
(6,1,1 



1,2,2 



16 

(4,2,1 
(4,1,2 



->■ 
->■ 

->■ 



8, 1) © (3, l) 2/3 © (3, l)_ 2 /3 © (1, l)o 
l,3)o 

1.1) 1 ©(1,1) ©(1,1)-1 

3, 2) 1/6 © (3, 2)_ 5/6 © (3, 2) 5/6 © (3, 2) 1/6 

3, l)_i /3 © (3, l)i /3 

1.2) 1/2 ©(l,2)_ 1/2 

3,2) 1/6 ©(l,2)_ 1/2 
3,1) 1/3 ©(3,1)_ 2/3 ©(1,1)i©(1,1)o 



(15) 



(16) 



(17) 



We have £7(1)?^ and C/(1)b_l from SU(2)r and 5'?7(4)c< respectively. Thus the weak 
hyper-charge generator (Ty) can be expressed as a linear combination of the generators of 
SU(2)r (T 3 n) and SU(A)c (Tb-l) sharing the same quantum numbers. In 10-dimensional 
representation T 3R , T B _ L and T Y are written as: 



T m = ^(0,0,0,0,0,0,^,-^,-^) 



(18) 



T^ = A^(-i, 44^, 0,0, 0,0) 



(19) 



/3 1 1 11111 11 1, 



Using these explicit forms of the generators we find the following relation 

T Y ~- 



^T m + \l T T B _ L 



I 5 V 5 

and this leads to the following mass relation, 



M 1 = - M 2R + - M 4C 



For 54-dimensional Higgs: 



(20) 



(21) 



(22) 



Using 54-dimensional Higgs we have [32], for D-parity even scenario, M^c = 1 and M 2 # 



M- 



2L 



). We have identified that M3 = M^c and M 2 = M 2 ^. Hence, using the above 



mass relation we obtain M 1 — — |. Therefore the gaugino mass ratio is given as: 

Mx : M 2 : M 3 = {-\) : {—) : 1 



(23) 
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We have already mentioned that the vev [13] of 770-dimensional Higgs can be expressed 
as a (45 x 45) diagonal matrix. So to calculate the gaugino masses, using 770-dimensional 
Higgs, we repeat our previous task in 45-dimensional representation. The explicit forms of 
Ty, T B _ L and T 3R are noted in the Appendix and find the same mass relation as in Equation 
(22) 

Mi = jj M 2R + H M iC (24) 

This mass relation among the U(l)y, SU(2) L and SU(3)c gauginos are independent of 
dimensions and representations. 

• 770-dimensional Higgs: 

Using 770-dimensional Higgs we find [13] for D-parity even case M 4 c = 2 and M 2R = M 2 l = 
5. Hence, using the above mass relation we obtain M\ = 3.8. Therefore the gaugino mass 
ratio is given as: 

Mi : M 2 : M 3 = 1.9 : 2.5 : 1 (25) 
We tabulate the gaugino mass ratios, obtained above, in Table 1. 



Representation 


M 3 : M 2 : M x at Mqut 


1 


1:1:1 


54: H -> 51/(4) x SU{2) x SU{2) 


l:(-3/2):(-l/2) 


770: H ->■ SU{A) x SU{2) x SC/(2) 


1:(2.5):(1.9) 



Table 1: i/ig/i sca/e gaugino mass ratios for the representations 54 and 770. 

Now we would like to address the question whether the D-Parity Odd case is phenomeno- 
logically viable. We know, if we use 210-dimensional Higgs instead of 54 or 770, the D-parity 
is broken by vev of 210 dimensional Higgs when 50(10) — > SU(A) C x SU{2) L x SU(2) R . 
This has a nice feature: M 2 l = —M 2R and M 4C = 0, which implies that at high scale M 3 
will be zero. We have analyzed the running of such high scale parameters numerically and 
found that these cases are not phenomenologically viable. With M 3 = at high scale, we 
will, always, hit Landau pole i.e. we will be in non-perturbative regime after running down 
to the low scale. This is not only the specialty of 210-dimensional Higgs, but also true for 
all the D -parity breaking scenarios. Hence we conclude that the D-parity non-conserving 
scenario doesn't fit into the non-universal gaugino mass framework, in particular, when we 
are bothered about a MSSM spectrum valid at the EWSB scale. 
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2.1 Implication of the Intermediate Scale 

We have an underlying assumption that the breaking of 5*0(10) GUT group to the inter- 
mediate gauge group and that to the SM takes place at the GUT scale itself, which is of 
course a simplification. But more interesting question to ask is how things will change if the 
intermediate scale is different from the GUT scale (which is usually the most realistic one)? 
Although we do not address this question in this analysis, we prefer to mention the crucial 
consequences of choosing an intermediate scale distinctly different from GUT scale: 

• In this case, the choice of the non-singlet Higgses will be restricted. Now, only those 
Higgses will contribute which have a singlet direction under the intermediate gauge 
group. Hence, the Higgses that break 50(10) directly to the SM at the GUT scale 
are disallowed. The choices of the non-singlet Higgses, in our analysis, are compatible 
even if the intermediate scale is different from the GUT scale. 

• The mass relation in Equation (22) is indeed independent of the intermediate scale as 
it is an outcome of purely group theoretical analysis. But the gaugino mass ratios will 
change depending on the choice of the intermediate scale due to the running of the 
gaugino masses from the GUT scale. 

3 Low energy spectra, Consistency and Benchmark 
Points 

Before we discuss in details the low-energy spectra for the non-universal inputs, we would like 
to mention a few points regarding the evolution of these gaugino mass ratios with different 
RGE specifications. As we know, in the one-loop RGE, the gaugino mass parameters do not 
involve the scalar masses [34], the ratios obtained at the low scale are independent of the high 
scale scalar mass input m . In addition, if we also assume no radiative corrections (R.C) 3 to 
the gaugino masses, the ratios at the EWSB scale are also independent of the choice of the 
gaugino mass parameters at the high scale. Instead, if one uses the two-loop RGE (scalars 
contributing to the gauginos), the values of the gaugino masses at the EWSB scale tend 
to decrease compared to the values obtained with one-loop RGE. Now, independently, the 
inclusion of R.C to the gaugino masses, makes the M 3 lower, but the values of M 1 and M 2 

3 The radiative corrections to the gauginos have been incorporated in calculating the physical masses after 
the RGE using the reference [35]. 
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become higher compared to the case of one- loop results with no R.C. When, one uses both 
the two-loop RGE and R.C to the gauginos, it is a competition between these two effects. 
In short, the gaugino mass ratios at the EWSB scale crucially depend on the choice of the 
RGE specifications. However, the dependence on the high scale mass parameters m and/or 
M 3 is very feeble 4 . 

We present in Table 2, the gaugino mass ratios at the EWSB scale for two different RGE 
conditions: 

• One- loop RGE with no R.C to the gaugino masses 



• Two- loop RGE with R.C to the gaugino masses 5 



Representation 


M-. : M 2 : M 1 
one- loop with No R.C 


M-. : M , : M 1 
two- loop with R.C 


1 (mSUGRA) 


1:0.27:0.13 


1:0.35:0.19 


54: H -> SU{A) x SU{2) x SU{2) 


l:(-0.40):(-0.06) 


l:(-0.55):(-0.10) 


770: H ->■ SU{A) x SU{2) x SU{2) 


1:0.67:0.24 


1:0.91:0.37 



Table 2: Low scale (EWSB) gaugino mass ratios for representations 54 and 770. 

The numerical results have been obtained using the spectrum generator SuSpect v2.3 
[36] with the pMSSM option. For the rest of our analysis we adhere to the second type of 
RGE specifications, two-loop RGE + R.C to the gauginos as mentioned earlier. The other 
broad specifications used for the scanning are listed below. 

• Full one-loop and the dominant two-loop corrections to the Higgs masses are incorpo- 
rated. 

• Gauge coupling constant unification at the high scale have been ensured and the corre- 
sponding scale has been chosen as the 'high scale' or 'GUT-scale' to start the running 
by RGE. All the non-universal inputs are provided at this scale using the pMSSM 
option. This is an artifact of choosing the intermediate scale set at the GUT scale 
itself. 

4 In particular, with change in M3 from 300-1000 GeV, with too=1000 GcV, the change in the ratios is 

within 10%, where as the ratios remain almost the same with change in m . 

5 The case with two-loop RGE + R.C to the gaugino masses have been obtained with mo = 500 

GcV. 
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• Electroweak symmetry breaking at the 'default scale' y/rn t ~ L m t ~ R has been set. 

• We have used the strong coupling a^(Mz) MS = 0.1172 for this calculation which is 
again the default option in SuSpect. 

• Throughout the analysis we have assumed the top quark mass to be 171.4 GeV. 

• All the scalar masses have been set to a universal value of mo and radiative electroweak 
symmetry breaking has been taken into account by setting high scale Higgs mass 
parameter M H 2 = M H 2 = m 2 and specifying sgn(fj,), which has been taken to be 
positive throughout the analysis. 

• All the trilinear couplings have been set to zero. 

• Tachyonic modes for sfermions and other inconsistencies in RGE, like Landau pole 
have been taken into account. 

• As we work in a -R-parity conserving scenario, stau-LSP regions have been identified 
as disfavoured. 

• Consistency with low-energy FCNC constraints such as those from b — > s 7 has been 
noted for each combination of the parameter space. We have used a 3a level constraint 
from b — > s 7 with the following limits [37]. 

2.77 x 1(T 4 < Br(b -> s 7 ) < 4.33 x 10~ 4 . (26) 

However, we must point out that we have taken all those regions as allowed where the 
value of b — > 57 is lower or within the constraint. 

• Regions allowed by all these constraints have been studied for the relic density con- 
straint of the cold dark matter (CDM) candidate (lightest neutralino in our case) and 
referred to the the WMAP data [15] within 3a limit 

0.091 < tt CDM h 2 < 0.128. (27) 

where flcDMh 2 is the dark matter relic density in units of the critical density and 
h = 0.71 ± 0.026 is the Hubble constant in units of 100 Km s" 1 Mpc" 1 . 

We have used the code microOMEGA v2.0.7 [38] for computing the relic density. 
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With these inputs, we scan the parameter space for a wide range of values of m and M 3 
6 for the non-universal gaugino mass ratios advocated above. 

The ratios obtained for 54 at the high scale (see Table 1) are actually the same as the one 
for the representation 24 in case of SU(5) (see [9, 10, 11, 24]). This observation differs from 
the earlier result available in [12]. The low-energy spectrum and its consistency for the case 
of 24 have been well-studied [21]. Without the inclusion of the intermediate breaking scale 
in case of S'O(IO), the case of 54 is difficult to distinguish from the one in SU(5). Anyway 
we do not address any such situations here and hence, refrain from illustrating the case of 
54. 

In Figure 1, we depict the results of the scan in the M 3 — itlq parameter space for the 
representation 770, i.e. breaking through G^2- Along the x-axis, high scale M 3 is varied 
from 100-2000 GeV and along the y-axis, high scale universal scalar mass itlq is varied in 
the same range. Our limit of the scan is motivated by the fact that we cover the low scale 
parameters well beyond the reach of the LHC The figure on the left hand side is for tan (3= 
5 and on the right hand side is for tan/3= 40. 

For tan/3= 5, full parameter space is allowed by REWSB, b — > 57 and other RGE 
constraints. The black region at the bottom (for M 3 = 400-1100 GeV and for very small 
values of m ) is disfavoured by the stau-LSP constraint. Hence, there is a large region of 
the parameter space (shown in red) which satisfies all the constraints and it is definitely 
within the reach of the LHC. We study the dark-matter constraints in this allowed region of 
parameter space and our conclusion is as follows: 

• For M 3 = 200 GeV, the allowed range of m spans around 200 GeV 

• For M 3 = 400 GeV, the allowed region is extremely narrow (because of the stau-LSP 
constraint) and is around m = 140 GeV 

• For M 3 = 600 GeV, m = 300 GeV is allowed 

• For M 3 = 800 GeV, the value of m goes as high as 1100 GeV 

We choose three benchmark points (BP1, BP2 and BP3, see Table 2 and 3) from here and 
study the collider signature. 

The figure on the right hand side of Figure 1 is with tan/3= 40 and is quite different 
from the one with tan/3= 5. In this figure, the red region is allowed by REWSB while the 
area in blue at small m is disfavoured by the stau-LSP constraint. Hence, here also, there 

6 Choice of M 3 automatically determines the values of Mi and M 2 for a choice of non-universality. 
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exists a large region of parameter space, sandwiched between these two, allowed by all the 
constraints for the study of dark-matter and collider search. Excepting for a very narrow 
region at the left bottom corner spanning 100-200 GeV of M3 or mo value (in green), the 
whole region is under the b — > sj upper limit. The dark matter study in this case, yields 
something special. We find almost all the regions to be lying below the lower bound of the 
WMAP data 7 . It is worthy to mention here that non-universal gaugino mass scenarios with 
M2 > M 3 at high-scale yield smaller value of /x generated from REWSB for a fixed M\ at 
TeV scale when compared with mSUGRA. This makes the neutralino-LSP more Higgsino 
like, increases the annihilation and consequently gives a larger region of the parameter space 
compatible with dark matter constraints [39]. We choose a couple of benchmark points (BP4 
and BP5, see Table 2 and 3) from here for the collider study. 




200 400 600 800 1000 1200 1400 1600 1800 2000 200 400 600 800 1000 1200 1400 1600 1800 2000 

M 3 M 3 



Figure 1: Parameter space scan for the representation 770 (breaking through G422) oi high 
scale M3 — mo plane (in GeV), showing the regions allowed by various constraints. Figure 
on the left hand side: tan/3= 5. The red region is allowed and the blue region at the 
bottom is disfavoured by stau-LSP. Figure on the right hand side: tan/3= Jf.0. The 
region in red is allowed by all constraints. Region in white is disfavoured by REWSB, the 
low-mo region, in blue, is disfavoured by stau-LSP and the left corner of the graph in green 
with small M 3 and m is disfavoured by b — >• 57. 



7 This is not strictly disfavoured as some other scenarios beyond the SM can co-exist and contribute. 
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Figure 2, shows similar parameter space scan for the case of mSUGRA, for tan/3= 5 
(left) and tan/3= 40 (right). These are presented to show the difference in the low-energy 
parameter space consistency patterns for different high scale gaugino mass inputs. These 
cases have been studied well and need not require much illustration. Our scan seems to match 
with the earlier available results [19] and show the robustness of our analysis. However, it is 
worth mentioning that the red region, excepting for the blue region at small m (disfavoured 
by stau-LSP), is allowed for tan/3= 5. For tan/3=40, a small region at the upper left 
corner with high values of m (1200-2000 GeV) and small M 3 (100-400 GeV) (in white) is 
disfavoured by REWSB. The blue region at the bottom is disfavoured for stau-LSP. 




200 400 600 800 1000 1200 1400 1600 1800 2000 200 400 600 800 1000 1200 1400 1600 1800 2000 



m 1/2 m 1/2 

Figure 2: Parameter space scan for mSUGRA at high scale m\ii — m plane (in GeV), 
showing the regions allowed by various constraints. Figure on the left side: tan/3 = 
5. The red region is allowed. The blue region at the bottom with low itlq is disfavoured by 
stau-LSP. Figure on the right side: tan/3= J^O. The red region is allowed. The region 
in white at the left corner with high values of mo (1200-2000 GeV) is disallowed by REWSB 
and the bottom of the figure in blue with low values of m (100-600 GeV), is disfavoured by 
stau-LSP. 

The benchmark points (BPs) chosen from Figure 1, to study the collider signature in 
context of the LHC, are presented in Table 3 and 4. In Table 3, we mention the high scale 
input parameters, while the low-energy spectra corresponding to these points have been 
mentioned in Table 4. The points have been chosen for two different values of tan /3, 5 (BP1, 
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BP2, BP3) and 40 (BP4, BP5) and have gluino masses around 500 GeV and 1000 GeV. 
These points for tan/3= 5 satisfy the WMAP data for the cold dark matter relic density- 
search, while the points for tan/3= 40 are all below the lower limit quoted by WMAP. The 
corresponding values of VtcDuh 2 have also been mentioned in Table 3. These points also 
obey the LEP bounds [40]. The model under scrutiny has been referred as 770-422 in Table 
3 and will be referred so in the following text. 



Benchmark Points 


Model 


M 3 


m 


tan (3 


^CDMh 2 


BP1 


770-422 


200 


200 


5 


0.124 


BP2 


770-422 


400 


140 


5 


0.125 


BP3 


770-422 


600 


300 


5 


0.126 


BP4 


770-422 


200 


200 


40 


0.0002 


BP5 


770-422 


400 


700 


40 


0.0162 



Table 3: Benchmark Points (BP): Models, High scale input parameters (in GeV), tan/3 

and VtcDuh 2 . 



Benchmark Points 


m g 


mq 12 


mr 

11,2 

m fl 




m xt 


n 


BP1 


499.95 


524 
305.26 


319.5 
246.79 


138.43 


203.5 
219.56 


221.4 


BP2 


938.8 


928 
551.83 


495 
313.93 


303.53 


375.96 
386.92 


381.18 


BP3 


1368.48 


1366 
815.14 


770.5 
513.84 


464.3 


515.65 
522.39 


515.27 


BP4 


499.5 


524.5 
315 


320 
172.25 


132.23 


170.45 
187.51 


178.36 


BP5 


966.58 


1145 
699.06 


857 
597.87 


246.73 


262.11 
269.58 


261.83 



Table 4 : Low-energy spectra for the chosen benchmark points (BP) (in GeV). 

In Table 4, we note the gluino mass (nig), average of the first two generation squark masses 
(rriq 12 ), average of the first two generation slepton masses (m ; ~ i2 ), lighter stau mass (m^), 
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lighter stop (m^), lightest neutralino (m^o), lighter chargino (m^±), 2nd lightest neutralino 
(m^o) as well as the value of fx, generated by REWSB at the BPs. It can also be noted that 
BP1, BP2, BP3 have bino dominated x?, while it is mixed and higgsino dominated in case 
of BP4 and BP5 respectively. The \\ * s mixed in case of BP1 and higgsino dominated in 
all the other cases. The lighter chargino is mostly higgsino dominated and degenerate with 
X21 excepting for the case of BP1. We would like to mention that, the composition (gaugino 
dominated, higgsino dominated or mixed) and the mass difference of the neutralinos and 
charginos get altered for different high scale gaugino non-universality. Given the similar 
values of squark and gluino masses in different non-universal schemes (actually similar choices 
of M3 and mo at the high scale), the electroweak gauginos become instrumental for a possible 
distinction between different GUT-breaking schemes, which might also get reflected at the 
collider signature in a favourable region of parameter space. 

To see such distinction in collider signature, we choose two points from mSUGRA scenario 
with suitable values of M1/2 such that the low scale gluino masses are around 500 GeV and 
1000 GeV (similar to the benchmark points (BP) selected above) in a region of parameter 
space that satisfies the WMAP data for cold dark matter constraint. It should be noted that 
once these points respect the CDM constraint, the value of m automatically get restricted 
for a particular choice of M 1 / 2 for a specific tan/3. Similar is the situation here where we 
have taken tan/3 = 5 for illustration. The chosen points are named as MSG1 and MSG2. 
The high scale parameters along with the VL C DMh 2 at these points are mentioned in Table 
5, while the corresponding low scale spectra are noted in Table 6. 



Points 


Model 


mi/2 


m 


tan /3 


&CDMh 2 


MSG1 


mSUGRA 


480 


100 


5 


0.111 


MSG2 


mSUGRA 


200 


70 


5 


0.128 



Table 5: mSUGRA points (MSG): Models, High scale input parameters (in GeV), tan/3 

and QcDMh 2 . 



Points 


m-g 


m qi,2 
™h 


m fl 


Xi 


m xf 

X2 


n 


MSG1 


1104.1 


992 
768.6 


272 
205.1 


195.66 


367.1 
367.4 


622.1 


MSG2 


493.25 


452 
323.6 


129 
101.3 


73.2 


131.4 
133.6 


281.4 
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Table 6 : Low-energy spectra for the chosen mSUGRA points (MSG). 

4 Collider Simulation and Numerical Results 

We would like to discuss the collider signature now, of the benchmark points advocated in 
the preceeding section. 

We first discuss the strategy for the simulation which includes the final state observables 
and the cuts employed therein. In the next subsection we discuss the numerical results 
obtained from this analysis. 

4.1 Strategy for Simulation 

The spectrum generated by SuSpect v2.3 as described in the earlier section, at the benchmark 
points are fed into the event generator Pythia 6.4.16 [41] by SLHA interface [42] for the 
simulation of pp collision with centre of mass energy 14 TeV. 

We have used CTEQ5L [43] parton distribution functions, the QCD renormalization and 
factorization scales being both set at the subprocess centre-of-mass energy \fs~. All possible 
SUSY processes and decay chains consistent with conserved .R-parity have been kept open. 
We have kept initial and final state radiations on. The effect of multiple interactions has been 
neglected. However, we take hadronization into account using the fragmentation functions 
inbuilt in Pythia. 

The final states studied here are : 

• Opposite sign dilepton (OSD) : (d ± (! T ) + (> 2) jets + Ejr 

• Same sign dilepton (SSD) : (£ ± £ ± ) + (> 2) jets + Ejr 

• Trilepton (3£ + jets): 3£ + (> 2) jets + Ejr 

• Hadronically quiet trilepton (31): 3£ + (0) jets + Ejr 

• Inclusive 4-lepton (4f): 4£ + X + Ejr 

where £ stands for final state isolated electrons and or muons, Ejr depicts the missing energy, 
X indicates any associated jet production. 

We will discuss these objects in details, that constitute the final state observables. The 
nomenclature assigned to the final state events in parantheses will be referred in the following 
text. 
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As defined in some earlier works [24], the absence of any jets with E T jet > 100 GeV 
qualifies the event as hadronically quiet. This avoids unnecessary removing of events along 
with jets originating from underlying events, pile up effects and ISR/FSR. The A£ events 
have been defined without putting an exclusive jet veto. 

Before we mention the selection cuts, we would like to discuss the resolution effects of the 
detectors, specifically of the ECAL, HCAL and that of the muon chamber, which have been 
incorporated in our analysis. This is particularly important for reconstructing Efc which is 
a key variable for discovering physics beyond the standard model. 

All the charged particles with transverse momentum, p? > 0.5 GeV 8 that are produced 
in a collider, are detected due to strong B-field within a pseudorapidity range \r]\ < 5, 
excepting for the muons where the range is \r]\ < 2.5, due to the characteristics of the muon 
chamber. Experimentally, the main 'physics objects' that are reconstructed in a collider, are 
categorised as follows: 

• Isolated leptons identified from electrons and muons 

• Hadronic Jets formed after identifying isolated leptons 

• Unclustered Energy made of calorimeter clusters with px > 0.5 GeV (ATLAS) and 
\r)\ < 5, not associated to any of the above types of high-.&r objects (jets or isolated 
leptons) . 

Below we discuss the 'physics objects' described above in details. 

• Isolated leptons (iso £): 

Isolated leptons are identified as electrons and muons with p T > 10 GeV and \rj\ <2.5. 
An isolated lepton should have lepton-lepton separation AR U >0.2, lepton-jet separation 
(jets with E T > 20 GeV) AR fj > 0.4, the energy deposit Y,E T due to \ow-E T hadron 
activity around a lepton within AR < 0.2 of the lepton axis should be < 10 GeV, where 
AR = \J Arj 2 + A0 2 is the separation in pseudo rapidity and azimuthal angle plane. The 
smearing functions of isolated electrons, photons and muons are described below. 

• Jets {jet): 

Jets are formed with all the final state particles after removing the isolated leptons from the 
list with PYCELL, an inbuilt cluster routine in Pythia. The detector is assumed to stretch 

8 This is specifically for ATLAS, while for CMS, p T > 1 GeV is used. 



19 



within the pseudorapidity range \r)\ from -5 to +5 and is segmented in 100 pseudorapidity 
(77) bins and 64 azimuthal ((f)) bins. The minimum E T of each cell is considered as 0.5 
GeV, while the minimum E T for a cell to act as a jet initiator is taken as 2 GeV. All the 
partons within Ai?=0.4 from the jet initiator cell is considered for the jet formation and 
the minimum J2 pa rton E T jet for a collected cell to be considered as a jet is taken to be 20 
GeV. We have used the smearing function and parameters for jets that are used in PYCELL 
in Pythia. 

• Unclustered Objects (Unc.O): 

Now, as has been mentioned earlier, all the other final state particles, which are not isolated 
leptons and separated from jets by AR >0.4 are considered as unclustered objects. This 
clearly means all the particles (electron/photon/muon) with 0.5 < E T < lOGeV and \rj\ < 5 
(for muon-like track \rj\ < 2.5) and jets with 0.5 < E T < 20GeV and \rj\ < 5, which are 
detected at the detector, are considered as unclustered energy and their resolution function 
have been considered separately and mentioned below. 

• Electron/Photon Energy Resolution : 



a(E)/E = a/\ r E®b®c/E 9 (28) 

Where, 

a = 0.03 [GeV 1 / 2 ], b = 0.005 & c = 0.2 [GeV] for \r)\ < 1.5 
= 0.055 = 0.005 = 0.6 for 1.5 < \rj\ < 5 



• Muon Pt Resolution : 



a(P T )/P T = a if P T < 100GeV (29) 

= a + b\og(P T /0 if P T > lOOGeV (30) 



Where, 

a= 0.008 & b= 0.037 for \r}\ < 1.5 
= 0.02 = 0.05 1.5 <M<2.5 



9 © indicates addition in quadrature 
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• Jet Energy Resolution 



a(E T )/E T = aj^E~ T (31) 

Where, 

a= 0.55 [GeV 1 / 2 ], default value used in PYCELL. 

• Unclustered Energy Resolution : 

a(E T ) = ay/EiEF™- * (32) 

Where, a ~ 0.55. One should keep in mind that the x and y component of E^ nc -° 
need to be smeared independently with same smearing parameter. 

All the smearing parameters that have been used are mostly in agreement with the 
ATLAS detector specifications and also have been discussed in details in [44] 

Once we have identified the 'physics objects' as described above, we sum vectorially 
the x and y components of the smeared momenta separately for isolated leptons, jets and 
unclustered objects in each event to form visible transverse momentum (pr)vis, 

(pr) ms = ^(E^) 2 + (Erf (33) 

where, J^Px = T,(Px)iso £ + T,(Px)jet + T,iPx)unc.o and similarly for J^Py We identify the 
negative of the (pt)vis as missing energy E{f. 

E{r=-{ PT ) ms (34) 

Finally the selection cuts that are used in our analysis are as follows: 

• Missing transverse energy Efr > 100 GeV. 

• Pt E > 20 GeV for all isolated leptons. 

• E T jet > 100 GeV and \r] jet \ < 2.5 

• For OSD, hadronically quiet trilepton (3£) and also for inclusive A£ events we have 
used, in addition, invariant mass cut on the same flavour opposite sign lepton pair as 
\M Z - M e+i _\ > 10 GeV. 
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We have checked the hard scattering cross-sections of various production processes with 
CalcHEP [45]. All the final states with jets at the parton level have been checked against 
the results available in [27]. The calculation of hadronically quiet trilepton rates have been 
checked against [46], in the appropriate limits. 

We have generated dominant SM events in Pythia for the same final states with same 
cuts, tt production gives the most serious backgrounds. We have multiplied the correspond- 
ing events in different channels by proper X-factor (= 2.23) to obtain the usually noted next 
to leading order (NLO) and next to leading log resummed (NLL) cross-section of tt produc- 
tion at the LHC, 908 pb (without taking the PDF and scale uncertainty), for m t around 
171 GeV [47]. The other sources of background include WZ production, ZZ production etc. 
The contribution of each of these processes to the various final states are mentioned in the 
Table 9. 

4.2 Numerical Results 

Figure 3 shows the effective mass distribution at the benchmark points and the corresponding 
mSUGRA ones in OSD events. Effective mass is defined as 

Effective mass = J2Miso t + £(Pr)jet« + E k (35) 

Figure 3 has been organised following the model inputs. Top left figure shows the distribu- 
tions at BP1, BP2 and BP3 chosen from 770-422 with tan/3=5, whereas the top right one 
contains BP4 and BP5 chosen from the same scenario with tan/3=40 and the one from tt 
production, the dominant process for the background. The bottom one is for mSUGRA, 
containing MSG1 and MSG2 with tan /3—5. The peak of the effective mass distribution 
corresponds to the threshold energy of the hard scattering process which is dominantly 
responsible for the final state under scrutiny. For OSD events, processes responsible are 
mostly the gg, gq and qq productions due to their SU(3) interactions, provided they are 
accessible to the LHC center of mass energy. In such cases, the threshold energy is around 
2m g or {m g + mq) or 2rriq. Now, in each of the BPs advocated here, m g = rriq and the thresh- 
old is approximately at 2m g . Our figures magnificently depict the correspondence with such 
threshold. For example, the peaks of BP1, BP4 and MSG2 are greater than 1000 GeV (where 
the gluino and squark masses are around 500 GeV). The reason that these distributions peak 
at higher values than the threshold, can be attributed to the fact that the final state con- 
sidered here, has a very large Ejp cut, (pr) cut on associated jets and leptons. While this 
indicates the robustness of our analysis, this also points to the deficiency in distinguishing 
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these non-universal models from the mSUGRA one with similar gluino masses. However, a 
possible way that could have been exploited is perhaps the effective mass distribution in 3£ 
events. This is expected as the dominant production process for this final state is X2 and 
Xi ± and these electroweak gauginos actually carry the information of different no n- universal 
gaugino mass inputs at the GUT scale. However, this was not very successful in our case 
due to small event rates. 

The missing energy distributions in OSD events at all the benchmark points have been 
shown in Figure 4. The organisation of the points remain the same as in Figure 3. In each 
case, the distribution starts from 100 GeV, as the event selection itself had this missing 
energy cut. As a result, all the points show a similar falling feature which indicates that the 
peak of the distribution is either small or around 100 GeV. The difference in the distributions 
is in the tail and is due to the hierarchy of the lightest neutralino masses. The heavier is 
the neutralino, the flatter is the distribution. Although this gives a nice distinction between 
the points with different gluino masses (and hence with different LSP masses), it is again, 
difficult to distinguish points with similar gluino masses. 

The numerical values of the event rates at the benchmark points are presented in Table 
7, while Table 8 contains the results in similar channels for the mSUGRA ones. We note 
the contributions to these channels from the SM background in Table 9. While we note 
that the results are widely different from each other for different BPs, we also point out 
the distinction with corresponding mSUGRA ones with similar gluino masses. For example, 
when we compare BP1, BP4 and MSG2 (all with gluino masses around 500 GeV), we note 
that the mSUGRA point yields much more events in almost all channels. This primarily 
has two reasons: one, the choice of the scalar mass parameter is very low for the mSUGRA 
one, compared to the non-universal case to obey CDM constraint and two, the non-universal 
scenario studied here, have higher values of Mi and M 2 at the high scale, which make the low- 
lying charginos and neutralinos heavier and correspondingly lower decay branching fraction 
through these to the leptonic final states. Similar observation can be made in an attempt to 
compare the event rates of BP2, BP3, BP5 and MSG1 {m~ g ~ 1000 GeV). The reason that 
BP2 has slightly higher event rates than MSG1 can be attributed to the fact that m g = 938.8 
GeV for BP1, which is smaller to be compared to MSG1 (m g = 1104 GeV). 

While we see that almost all the channels at the BPs rise sufficiently over the background 
fluctuations, the hadronically quiet trilepton channel gets submerged in to the background 
even for an integrated luminosity 100/6 -1 at the points BP2, BP3, BP5 and MSG1. This 
is because of the very high values of m^ 2 o and m^ x ± with very small Am( x ~ 2 o/^ 1 ±_ x ~ 1 0). The 



23 




Effective Mass (GeV) 

Figure 3: Effective mass distribution in OSD events. Top left: BP1 (Red), BP2 (Blue) 
and BPS (Pink,) chosen from 770-422 with tan/3=5 have been plotted. Top right: BP 4 
(Red) and BP 5 (Blue) chosen from 770-422 with tan f3 =40 and ti (T'inkj have been plotted. 
Bottom: MSG1 (Red) and MSG2 (Blue) chosen from mSUGRA with tan/3=5 have been 
plotted. CTEQ5L pdfset was used. Factorisation and Renormalisation scale has been set to 
t^F = t^R = V§, sub-process centre of mass energy. 
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Missing Energy (GeV) 

Figure 4: Missing energy distribution in OSD events. Top left: BP1 (Red), BP2 (Blue) 
and BPS (Pink,) chosen from 770-422 with tan/3=5 have been plotted. Top right: BP 4 
(Red) and BP 5 (Blue) chosen from 770-422 with tan j3 =40 and ti (Pink) have been plotted. 
Bottom: MSG1 (Red) and MSG2 (Blue) chosen from mSUGRA with tan/3=5 have been 
plotted. CTEQ5L pdfset was used. Factorisation and Renormalisation scale has been set to 
— — V§, sub-process centre of mass energy. 
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significance in most of these points, in most of the channels (excepting the 3£), are so high 
that it is very unlikely to be affected by the systematic errors. We would also like to point out 
that all the channels in BP1, BP4 and MSG2 rise over the background even for an integrated 
luminosity of 10/6 -1 (only exception being the a S e for BP4), while others are suppressed by 
the background. In Table 10, we summarise this information for each of the channels and 
parameter points, for an integrated luminosity of 

In the tables, the cross-sections are named as follows: cr OSD for OSD, cfssd 

for SSD, 

&u+jets f° r (3£ + jets), (T31 for (3£) and for inclusive 4 lepton events 4£. 



Benchmark Points 


&OSD 


0~SSD 


0~3t+jets 


0"3£ 


0~4£ 


BP1 


597.7 


102.9 


57.6 


13.10 


10.30 


BP2 


70.1 


15.7 


11.6 


0.09 


1.64 


BP3 


17.9 


2.2 


4.0 


0.01 


1.16 


BP4 


398.4 


120.8 


30.8 


4.50 


3.28 


BP5 


26.5 


12.2 


3.9 


0.05 


0.33 



Table 7: Event-rates (fb) in multilepton channels at the chosen benchmark points. CTEQ5L 
pdfset was used. Factorisation and Renormalisation scale has been set to hf = Hr = 
subprocess centre of mass energy. 



mSUGRA Points 


&OSD 


0~SSD 


&M+jets 




&u 


MSG1 


50.9 


16.4 


10.6 


0.08 


1.34 


MSG2 


2781.8 


175.9 


285.1 


12.26 


99.89 



Table 8: Event-rates (fb) in multilepton channels at the mSUGRA points. CTEQ5L pdfset was 
used. Factorisation and Renormalisation scale has been set to fi F = /j, R = y/§, subprocess 
centre of mass energy. 



SM Processes 


0~OSD 


0~SSD 


°~3£+jets 


0~3£ 




tt 


1102 


18.1 


2.7 


5.3 


0.0 


ZZ,WZ,ZH,Z-f 


16.3 


0.3 


0.5 


1.1 


0.4 


Total 


1118.3 


18.4 


3.2 


6.4 


0.4 



Table 9: Event-rates after cut (fb) in multilepton channels from the dominant SM back- 
grounds. The event rates in different channels for tt production have been multiplied by by 
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proper K -factor (2.23) to obtain the usually noted NLO+NLL cross-section oftt [47]. CTEQ5L 
pdfset was used. Factorisation and Renormalisation scale has been set to pip — — 
subprocess centre of mass energy. 



Model Points 


OSD 


SSD 


3£ + jets 


u 


4£ 


BP1 


V 


V 




V 




BP2 






V 


X 




BP3 


X 


X 


V 


X 




BP4 


V 


V 


V 


V 




BP5 


V 






X 


X 


MSG1 


V 


V 


V 


X 


V 


MSG2 


V 




V 


V 





Table 10 : 5a visibility of various signals for an integrated luminosity of 30 fb 1 . A 
indicates a positive conclusion while a x indicates a negative one. 

We also compare these results in the ratio space of events which is demonstrated in Figure 
5 in form of barplot. The advantage of going to the ratio space is the uncertainties due to 
the choice of pdfsets, jet energy scale get reduced. Here we take the ratios of all events with 
respect to the OSD and referred as SSD/OSD, 3L+JETS/OSD, 3L/OSD and 4L/OSD along 
the x-axis of the barplot. As earlier, we divide the BPs and the MSGs in two categories: one, 
with BP1, BP4 and MSG2 (trig ~ 500 GeV), which is shown in the left panel of the Figure 
5; two, BP2, BP3, BP5 and MSG1 [m g ~ 1000 GeV), which is shown in the right panel of 
the Figure 5. We note that BP1, BP4 and MSG2 are well distinguished from SSD/OSD, 
3L/OSD and 4L/OSD. While BP2 and MSG1 can not be distinguished very well from each 
other, identification of BP3, BP5 and MSG1 is quite apparent from SSD/OSD and 4L/OSD 
events. 

In a nutshell we can summarise that, it is indeed possible to distinguish the non-universal 
gaugino mass scenario advocated here, from the mSUGRA ones with similar gluino masses. 
This is possible in both the absolute event rates or from the ratio plots shown here. However, 
the distinguishability reduces with increasing gluino mass. 
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Event Ratios Event Ratios 

Figure 5: Event ratios of different benchmark points compared with mSUGRA. Along the x- 
axis, ratios of events with respect to OSD have been taken. The relative values of these ratios 
for different benchmark points are indicated along the y-axis. Figure on the left: MSG2 
(Red), BP1 (Gieen) and BP4 (Blue) are compared, nig is around 500 GeV in each case. 
3L/OSD and 4L/OSD have been multiplied by 10. Figure on the right: MSG1 (Red), 
BP2 (Green), BPS (Blue ) and BPS (Pinkj are compared, nig is around 1000 GeV in each 
case. 3L/OSD has been multiplied by 100 and 4L/OSD has been multiplied by 10 to accom- 
modate them in the same figure. CTEQ5L pdfset was used. Factorisation and Renormalisation 
scale has been set to u, F — u, R — vl, sub-process centre of mass energy. 
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5 Summary and Conclusions 



We have derived non-universal gaugino mass ratios for the representations 54 and 770 for the 
breaking chain SU{A) C x SU(2) L x SU(2) R (O422) in a 50(10) SUSY-GUT scenario. We have 
assumed that the breaking of 50(10) to the intermediate gauge group and the latter in turn 
to the SM gauge group takes place at the GUT scale itself. We point out some errors in the 
earlier calculation and derive new results on the gaugino mass ratios. We scan the parameter 
space with different low-energy constraints taken into account and point out the allowed 
region of the parameter space. We also study the dark matter constraint in these models 
and study collider simulation at some selected benchmark points in context of the LHC. The 
scans presented, have many interesting features that might help us in understanding the 
correlation between high scale input and low-energy spectra. We must mention here that 
the study is limited by the assumption that the series of symmetry breaking occurs at the 
GUT scale itself. It is essentially a simplification, although we know that the mass relation 
in Equation (22) is not going to change with different choices of the intermediate scale, 
while the running of the masses from the GUT scale to the intermediate scale will eventually 
change the gaugino mass ratios. However, our choice of the non-singlet Higgses, apart from 
conserving Z)-parity, are compatible even if the intermediate scale is different from the GUT 
scale. Within such a framework we have performed a collider study which is more illustrative 
than exhaustive. It nonetheless elicits some characteristics of the signature space for such 
high scale ratios in context of the LHC. It is worth mentioning in this context that the 
comparison between the non-universal inputs with the mSUGRA ones yields a significant 
distinction in the multilepton channel parameter space in context of the LHC. This feature 
might be important in pointing out the departure of the 'signature space' with the inclusion 
of non-universality at the high scale. 

While this paper was in preparation, we came across the reference [48], where the issue 
of gaugino mass non-universality in the context of 50(10) and E(G) has been addressed. 
While we agree completely with the corrected gaugino mass ratios for G 42 2, our analysis 
in addition, point out the phenomenologically viable situations where the choices of the 
Higgses get restricted with D-parity conservation and inclusion of the intermediate scale 
different from the GUT scale. Furthermore, the low-energy particle spectra in different 
cases have been derived in a comparative manner and the allowed regions of the parameter 
space consistent with low-energy and dark matter constraints are obtained in each case. In 
addition, we have predicted event-rates for such a breaking chain in a multichannel study 
pertinent to the LHC. The distinguishability of relative rates in different channels has also 
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been explicitly demonstrated by us. 
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Appendix 



Using Equation (15) in Section 2, T Y , T B -l and T 3R in 45-dimension can be written as: 

T Y = \^diag(0, ..0, -, -, -, --, --, --, 0, ..0, 1, 0, -1, -, .., -, --, .., --, -, .., -, --, .., --)(36) 
Y V5 yv ^^'3'3'3' 3' 3' ' ' '6' J 6 \ 6' J 6 ' 6' J_6 \ 6 J 6, ; 



3, , 444444 2 
7b_l = W -diagQJQ, -, -, -, --, --, --, 0^0, --, 
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(37) 



T 3i? = diagi^^O, 1, 0, -1, 1 .., ± ~, .., -± ± .., 1 ~, .., -±) (38) 

15 3 6 6 6 6 

These lead to the following relation 

T Y = J^T 3R + J^T B _ L (39) 
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